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^ . Abstract 



We discuss charge symmetry and charge independence breaking in an effective 
field theory approach for few-nucleon systems. We systematically introduce strong 
isospin-violating and electromagnetic operators in the theory. The charge depen- 
C^ \ dence observed in the nucleon-nucleon scattering lengths is due to one-pion ex- 

change and one electromagnetic four-nucleon contact term. This gives a parame- 
ter free expression for the charge dependence of the corresponding effective ranges, 
which is in agreement with the rather small and uncertain empirical determinations. 
We also compare the low energy phase shifts of the nn and the np system. 
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1. It is well established that the nucleon-nucleon interactions are charge dependent (for 
a review, see e.g.|]I|]). For example, in the ^5*0 channel one has for the scattering lengths 
a and the effective ranges r {n and p refers to the neutron and the proton, in order) 

'^«CIB = 2 ('^"" + ""PP^ ~ ^^P ^ ^-^ ^ ^-^ ^"^ ' 

ArciB = - (r„„ + rpp) -r„p = 0.05 ± 0.08 fm. (1) 

These numbers for charge independence breaking (CIB) are based on the Nijmegen po- 
tential and the Coulomb effect for pp scattering is subtracted based on standard methods. 
The charge independence breaking in the scattering lengths is large, of the order of 25%, 
since a„p = (—23.714 ±0.013) fm. In addition, there are charge symmetry breaking (CSB) 
effects leading to different values for the pp and nn threshold parameters, 

AacsB = a-pp - ann = 1-5 ± 0.5 fm , 

ArcsB = rpp - r„„ = 0.10 ± 0.12 fm. (2) 

Both the CIB and CSB effects have been studied intensively within potential models of the 
nucleon-nucleon (NN) interactions. In such approaches, the dominant CIB comes from 
the charged to neutral pion mass difference in the one-pion exchange (OPE), Aa^f^ ~ 
3.6 ± 0.2 fm. Additional contributions come from 77? and 2tt (TPE) exchanges. Note also 
that the charge dependence in the pion-nucleon coupling constants in OPE and TPE 
almost entirely cancel. In the meson-exchange picture, CSB originates mostly from p — u 
mixing, Aa^gg ~ 1.2 ± 0.4 fm. Other contributions due to n — rj, n — t]' mixing or the 
proton-neutron mass difference are known to be much smaller. 

Within QCD, CSB and CIB are of course due to the different masses and charges of 
the up and down quarks. Such isospin violating effects can be systematically analyzed 
within the framework of chiral effective field theories. In the two-nucleon sector, a com- 
plication arises due to the unnaturally large S-wave scattering lengths. This can be dealt 
with in two ways. One is the "hybrid" approach due to Weinberg 0, in which chiral 
perturbation theory is applied to the interaction kernel sewed to realistic nuclear wave 
functions obtained by conventional means. This approach has been successfully applied to 
a variety of processes which are dominated by OPE, a particularly striking one being the 
prediction of the electric dipole amplitude for neutral pion production off deuterons |§] 
recently measured at SAL [Q. A different and more systematic fashion to deal with the 
unnaturally large scattering lengths is the recently proposed power divergence subtrac- 
tion scheme (PDS) proposed by Kaplan, Savage and Wise (KSW) 0.|3 Essentially, one 
resums the lowest order local four-nucleon contact terms ~ Co(A^^A^)^ (in the S-waves) 
to generate the large scattering lengths and treats the remaining effects perturbatively, in 
particular also pion exchange. This means that most low-energy observables are domi- 
nated by contact interactions. The chiral expansion for NN scattering entails a new scale 
AjVAT of the order of 300 MeV, so that one can systematically treat external momenta 
up to the size of the pion mass. There have been suggestions that the radius of conver- 
gence can be somewhat enlarged [^], but in any case A^^ is considerably smaller than 

"^■^There exist by now modifications of this approacli and it as been argued tliat it is equivalent to 
cut-ofF schemes. We do not want to enter this discussion here but rather stick to its original version. 



the typical scale of about 1 GeV appearing in the pion-nucleon sector. A status report of 
the various observables calculated within this framework can be found in ref.0. In this 
context, it appears to be particularly interesting to study CIB (or in general isospin vio- 
lation) which is believed to be dominated by long range pion effects. That is done here 



m 



First, we write down the leading strong and electromagnetic four-nucleon contact terms. 
It is important to note that in contrast to the pion or pion-nucleon sector, one can not 
easily lump the expansion in small momenta and the electromagnetic coupling into one 
expansion but rather has to treat them separately. Then we consider in detail CIB. The 
leading effect starts out at order aQ'"^, where Q is the generic expansion parameter in the 
KSW approach. It stems from OPE plus a contact term of order a with a coefficient of 
natural size that scales as Q"^- Similarly, the leading CSB effect are four-nucleon contact 
terms of order a and order rriu — niii, which also scale as Q^^. While in the case of Eq 
this scaling property is enforced by a cancellation of a divergence, the situation is a priori 
different for CSB. However, for a consistent counting of all isospin breaking effects related 
to strong or em insertions, one should count the quark mass difference and virtual photon 
effects similarly. Note, however, that these CIB and CSB terms are numerically much 
smaller than the leading strong contributions which scale as Q~^ because a -C 1 and 
{niu — m^)/ A^ <^ 1. The corresponding constants, which we call Eq , together with the 
strong parameters (as given in the work of KSW) can be determined by fitting the three 
scattering lengths app, ann, cinp and the np effective range .p^ That allows to predict the 



momentum dependence of the np and the nn ^Sq phase shifts. Based on these observation, 
we can in addition give a general classification for the relevant operators contributing to 
CIB and CSB in this scheme. Additional work related to long-range Coulomb photon 
exchange is necessary in the proton-proton system. We do not deal with this issue here 



but refer to recent work using EFT approaches in refs. [|TT| ||T2 



2. First, we discuss the various parts of the effective Lagrangian underlying the analysis 
of isospin violation in the two-nucleon system. To include virtual photons in the pion 
and the pion-nucleon system is by now a standard procedure |[T^-[|l^. The lowest order 
(dimension two) pion Lagrangian takes the form 

/:.. = ^(v^f/v^^yt + xu^ + x^u) + c{ququ^) , (3) 

with /,r = 92.4 MeV the pion decay constant, V^ the (pion) covariant derivative containing 
the virtual photons, ( ) denotes the trace in flavor space, x contains the light quark masses 
and the last term, which contains the nucleon charge matrix (5=ediag(l, O),0 leads to 
the charged to neutral pion mass difference, 5m^ = m^± — wi^o, via 6m'^ = SnaC/f^, i.e. 
C = 5.9 ■ 10~^ GeV^. Note that to this order the quark mass difference rriu — rrid does 
not appear in the meson Lagrangian (due to G-parity). That is chiefly the reason why 
the pion mass difference is almost entirely an electromagnetic (em) effect. The equivalent 

'^''For a first look at these effects in an EFT framework, see the work of van Kolck [g[||^. Electromagnetic 
corrections to the one-pion exchange potential have been considered in |10|| . 
^^Whenever we talk of the pp system, we assume that the Coulomb effects have been subtracted. 
^^Or equivalently, one can use the quark charge matrix e(l + tj^JI. 



pion-nucleon Lagrangian to second order takes the form 

^ 4 V^ ^ am) ^^'' ^j]'"*^J' + ^5 [X+ - 2^X+)j +---JN , (4) 

which is the standard heavy baryon effective Lagrangian in the rest-frame w^ = (1, 0, 0, 0). 
M is the nucleon mass and m^ the chiral viel-bein, m^ ~ —idf^(f)/ fj^ + . . . . The four- 
nucleon interactions to be discussed below do not modify the form of this Lagrangian (for 



a general discussion, see e.g. ref|]20[). Strong isospin breaking is due to the operator ~ C5. 



Electromagnetic terms to second order are given by |T8| 

CTn = flN^ {h{Ql - Q'-) + /2g+(g+) + UQl + Ql) + /4(Q+)'} n , (5) 

with Q± = uQ'^u ± u^Qu^ and A = A — {A)/2 projects onto the off-diagonal elements of 
the operator A. Evidently, the charge matrices always have to appear quadratic since a 
virtual photon can never leave a diagram. The last two terms in eq. (^) are not observable 
since they lead to an equal em mass shift for the proton and the neutron, whereas the 
operator ~ /2 to this order gives the em proton-neutron mass difference. In what follows, 
we will refrain from writing down such types of operators which only lead to an overall 
shift of masses or coupling constants. We note that in the pion and pion-nucleon sector, 
one can effectively count the electric charge as a small momentum or meson mass. This is 
based on the observation that M^/A^ ~ e/y/Air = ^/a ~ 1/10 since A^ ~ 47r/^ = 1.2 GeV. 
It is thus possible to turn the dual expansion in small momenta/meson masses on one side 
and in the electric coupling e on the other side into an expansion with one generic small 
parameter. We also remark that from here we use the fine structure constant a = e'^/Air 
as the em expansion parameter. 

We now turn to the two-nucleon sector, i.e. the four-fermion contact interactions 
without pion fields. Consider first the strong terms. Up to one derivative, the effective 
Lagrangian takes the form 

Cfj, = h{N^Ny + k{N^aNY + hiNHx+)N)iN^N) +hiN^X+N)iN^N) 

+ kiN^a{x+)N){N^aN) + k{N^ax+N){N^aN) + . . . , (6) 

where the ellipsis denotes terms with two (or more) derivatives acting on the nucleon 
fields. Similarly, one can construct the em terms. The ones without derivatives on the 
nucleon fields read 

C^j, = N^{n{Ql-Q'_) + r,Q4Q^)}N{N^N) 

+ N^a{r,{Ql-Ql)+nQ+{Q+)}N{N^aN) 

+ Art {r,Q^ + r,{Q+)} N{N^Q+N) + N^a {rjQ+ + rs{Q+)} N{N^aQ+N) 

+ r,iN^Q+N)'' + noiN^aQ+N)'' . (7) 

There are also various terms resulting form the insertion of the Pauli isospin matrices r in 
different N'^N binomials. Some of these can be eliminated by Fierz reordering, while the 



others are of no importance for our considerations. Note that from now on we consider 
em effects. The leading CSB ~ my_ — ma has the same structure as the corresponding em 

(2) 

term and thus its contribution can be effectively absorbed in the value oi Eq , as defined 
below. We remark since A^vat is significantly smaller than A^, it does not pay to treat 
the expansion in the generic KSW momentum Q simultaneously with the one in the fine 
structure constant (as it is done e.g. in the pion-nucleon sector). Instead, one has to assign 
to each term a double expansion parameter Q^a™, with n and m integers. Lowest order 
charge independence breaking is due to a term ~ (NW^N)"^ whereas charge symmetry 
breaking at that order is given by a structure ~ {NW^N){N^N). In the KSW approach, it 
is customary to project the Lagrangian terms on the pertient NN partial waves. Denoting 
by /? the ^5*0 partial wave for a given cms energy Ecms, the Born amplitudes for the lowest 
order CIB and CSB operators between the various two-nucleon states takes the form 



{(3,pp\L^;rN\(^,pp) = -(^^^a[E^^^ + Ef^ 



{(3,nn\C^^\(3,nn) = - (^Ml^ a {e^,'^ - E^ 



(1 21 

where we will determine the coupling constant E^^ ' later on and also derive the scaling 
properties of the £^2n ■ The terms with the superscript '(!)' refer to CIB whereas the 
second ones relevant for CSB are denoted by the superscript '(2)'- Higher order operators 
are denoted accordingly. There is, of course, also a CIB contribution to the np matrix 
element. To be consistent with the charge symmetric calculation of ref. 0], we absorb its 
effect in the constant D2, i.e. it amounts to a finite renormalization of D2 and is thus not 
observable. In eq.(^), p = yJMEcms is the nucleon cms momentum. 

3. Consider now the effect of the charged to neutral pion mass difference 5m? , see the 
upper left diagramm in fig. 1. OPE between two neutrons or two protons can obviously 
involve charged and neutral pions. The mass difference can be treated in two ways. As 



proposed in ref.|T6[, one can modify the pion propagator, 

^ ^ ^ " [P- ml, -5m'' [1-5^-)] ' " ~lf~ ' ^' 

with i the pion four-momentum and 'a, h' isospin indices. Since we are interested only in 
the leading corrections ~ 5m? ~ a, it suffices to work with the expanded form of eq.(P), 

From this we conclude that OPE diagrams with different pion masses have the isospin 
structure O12 = r(i)aA"*r(i){, = (a + 6) r(i) f(2) — ^'^m '''{2) ^^^ ^^^^ to CIB since 

{pp\0i2\pp)cs = {nn\0i2\nn) = - , 
{np\0i2\np) + {np\0i2\pn) = - + - , (11) 



for the various isospin components of the two-nucleon system and 'Cs' stands for Coulomb- 
subtracted. Obviously, these effects are of order aQ~^. The np amplitude was already 
calculated by KSW. We have worked out the leading corrections A^ = Ann — ^np = 
App — Anp due to the pion mass difference. The pertinent diagrams are shown in fig. 1. 
We follow the notation of KSW in that we call these corresponding three amplitudes 
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where the first (second) subscript refers to the power in a (Q) and the super- 



scripts to the first three diagrams of the figure. We find 
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(12) 



with A-\ = Ao-i the leading term in the expansion of the np ^Sq amplitude 



A-i 



Co(/i) 



l + Co(/i)M(/i + ip)/47r 



(13) 



Here, /x is the PDS regularization scale. Note while the diagrams II and III are convergent, 
IV diverges logarithmically. Therefore, the Lagrangian must contain a counterterm of 
the structure Eq {ix)a{NW^N){NW^N) (cf. section 2) since it is needed to make the 
amplitude scale-independent. Note that for graph IV we have used the same subtraction 
as performed in [Q . Consequently, for operators of this type with 2n derivatives we can 
establish the scaling property £'2^ ~ Q"^^". This does not contradict the KSW power 
counting for the isospin symmetric theory since a <^ 1. Stated differently, the leading CIB 
term of order aQ""^ is numerically much smaller than the strong leading order contribution 
~ Q~^. The insertion from this contact term is shown in the last diagram of fig. 1 and 
leads to an additional contribution to A^. In complete analogy, we can treat the leading 
order CSB effect which is due to an operator of the form aE^ {NW3N){N'^N). This term 
is, however, finite. Putting pieces together, we get 
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K^ + 4 
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where the coupling constants Eq ' (yu) obeys the renormalization group equations. 



(14) 
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(15) 



Note that from here on we do no longer exhibit the scale dependence of the various 

(1 2) 

couplings constants Eq ' ,Co,2,-D2- We can now relate the pp and nn scattering lengths 
to the np one (of course, in the pp system Coulomb subtraction is assumed), 
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+ A, 



+ A, (16) 



^g\{~CoM{m - 2^) + Stt + CoMmln(mV/i2)) 
167rmCo/2 

For the effective ranges, we have only CIB 

_ _ , .2 9\{CoM^, + 47r)(CoM(3/x - 2m) + 127r) 
rnn-Tpp-rnp + Om QT^Mm^Clfl ' 

Note that this last relation is scale-independent and that it does not contain any unknown 
parameter. We remark that for the CIB scattering lengths difference the pion contribution 
alone is not scale-independent and can thus never be uniquely disentangled from the 
contact term contribution r^ Eq . While the leading OPE contribution resembles the 
result obtained in meson exchange models, the mandatory appearance of this contact 
term is a distinctively new feature of the effective field theory approach. It is easy to 
classify the leading and next-to-leading em corrections to these results. At order aQ~^, 
one has the contribution from two potential pions with the pion mass difference and 
also contact interactions with two derivatives. Effects due to the charge dependence of 
the pion-nucleon coupling constants, i.e. isospin breaking terms from £^™ , only start 
to contribute at order aQ^. Such effects are therefore suppressed by two orders of Q 
compared to the leading terms. This finding is in agreement with the various numerical 
analyses performed in potential models. We now turn to CSB. Here, to leading order 
there is simply a four-nucleon contact term proportional to the constant Eq . Its value 
can be determined from a fit to the empirical value given in eq.(|^). First corrections to 
the leading order CSB effect are classified below. 

4. We now turn to the numerical analysis considering exclusively the ^5*0 channel. At 
leading order, all parameters can be fixed from the pertinent scattering length. Already 
in ref. it was pointed out that there are various possibilities of fixing the next-to- 
leading order constants. One is to stay at low energies (i.e. fitting a and r) or performing 
an overall fit to the phase shift up to momenta of about 200 MeV. In this latter case, 
however, the resulting low energy parameters deviate from their empirical values and also, 
at p ~ 150 MeV, the correction is as big as the leading term. Since we are interested in 
small effects like CIB and CSB, we stick to the former approach and use the ^5*0 np phase 
shift around p = to fix the parameters as shown by the dashed curve in fig. 2. More 
precisely, we fit the parameters Co,2, -D2 to the np phase shift under the condition that the 
scattering length and effective range are exactly reproduced. The two new parameters 

(1 2) 

Eq ' are determined from the nn and pp scattering lengths. The resulting parameters at 



H = m are of natural size, 

Co = -3.46 fm^ , C2 = 2.75 fm^ , D2 = 0.07 fm^ , 
eJ^) = -6.47 fm^ , E^^^ = 1.10 fm^ . (18) 

To arrive at the curves shown in fig. 2, we have used the physical masses for the proton 

(1 2) 
and the neutron. We stress again that the effect of the em terms ~ -Eg is small because 

of the explicit factor of a not shown in eg . (pISD . Having fixed these parameters, we can 

now predict the nn ^So phase shift as depicted by the solid line in fig.2. It agrees with 

nn phase shift extracted from the Argonne V18 potential (with the scattering length and 

effective range exactly reproduced) up to momenta of about 100 MeV. The analogous 

curve for the pp system (not shown in the figure) is close to the solid line since the CSB 

effects are very small. 

5. Finally, we can give a summary of the various leading (LO), next-to-leading (NLO) 
and next-to-next-to-leading order (NNLO) contributions to CIB and CSB, with respect 
to the expansion in Q, to leading order in a and the light quark mass difference. Consider 
first CIB. 

LO, aQ^^ Electromagnetic corrections to one-pion exchange (pion mass difference) and 
electromagnetic four-nucleon contact interactions with no derivatives. 

NLO, aQ"^, eQ~^ Em corrections to two-pion exchange, four-nucleon contact terms with 
two derivatives and insertions proportional to the strong neutron-proton mass split- 
ting. 

NNLO, aQ^ Em corrections to three-pion exchange and four-nucleon contact terms with 
four derivatives as well as charge dependent coupling constants from the electromag- 
netic pion-nucleon Lagrangian. 

Note that there are no CIB effects due to the light quark mass difference linear in e = 
TUu — m^pjl except the trivial effects related to the strong pn mass difference. We now 



turn to CSB. The pattern of the various contributions looks different: 

LO, aQ~'^,eQ~'^ Electromagnetic and strong isospin-breaking four-nucleon contact in- 
teractions with no derivatives. 

NLO, aQ^^,eQ^^ Em four-nucleon contact terms with two derivatives, strong isospin- 
breaking contact interaction with two derivatives and insertions proportional to the 
neutron-proton mass difference. 

Here, we assume the same scaling in powers of Q for the em and the strong coupling 
constants. We remark that the leading order CSB effects do not modify the effective 
range but the corresponding scattering length. 



"^^One could also use the dimensionsless small parameter e = (tou — md)/{mu + md) — 1/3, see e.g. 
ref. |B|. We prefer to keep the notation which is most commonly used in the pion-nucleon sector, see e.g. 
ref.|2l|. 



6. In summary, we have considered electromagnetic and strong isospin violation in low- 
energy nucleon-nucleon scattering in the effective field theory formalism developed in 
ref.0. In particular, the leading charge independence breaking effect is due to a com- 
bination of the neutral to charged pion mass difference in one-pion exchange diagrams 
together with an electromagnetic four-nucleon contact term. Its corresponding coupling 
constant scales as Q~^ but is numerically suppressed by the explicit appearance of the 
fine structure constant a ~ 1/137. We have shown how the KSW power counting has to 
be modified in the presence of isospin violating operators. We explicitely evaluated the 
^5*0 phase shifts for the np, nn and Coulomb-subtracted pp systems at next-to-leading 
order. In addition, we have given a general classification of the various CIB and CSB 
corrections. It would be interesting to extend this formalism to other partial waves and 
to higher energies so as to investigate e.g. isospin violation in pion production. 
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Figure 1: Relevant graphs contributing to charge independence breaking at leading order aQ~^. 
The blob stands for the resummation of the lowest order {N^ N)"^ contact terms ~ Cq. The open 
(filled) circle denotes a pion mass insertion ~ dm? (an insertion of the leading four-nucleon 
operators ~ oEq ). For charge symmetry breaking, the leading contribution is given by the last 

(2) 

diagram with the filled circle denoting an insertion~ ctEk . 
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Figure 2: ^5*0 phase shifts for the np (dashed line) and nn (solid line) systems versus the 
nucleon cms momentum. The emipirical values for the np case (open squares) are taken from 



the Nijmegen analysis [22|. The open octagons are the nn "data" based on the Argonne V18 
potential. 
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